clear all
clc
close all

Rz=6400e3;
Rorb=19100e3;
c=3e8;
T=0:60:24*3600;
w=2*pi/(24*3600);
Rr=[0;Rz];
N=size(T,2);
Rsv=[255e5;0];
Rv=[255e5;-64e5];
f0=1602e6;

for i=2:N
   
Rsv(:,i)=[(Rz+Rorb)*cos(w*T(i)); (Rz+Rorb)*sin(w*T(i))];
Rv(:,i)= Rsv(:,i)-Rr;
Vsu(:,i)=(Rv(:,i)-Rv(:,i-1))/60;
Vu(:,i)=-dot(Vsu(:,i),Rv(:,i))/norm(Rv(:,i));
fd(:,i)=f0*Vu(:,i)/(c);
end;



plot(Rsv(1,:),Rsv(2,:),Rv(1,:),Rv(2,:))
plot(T,fd)



